An Efficient Approach for Multifrontal Al- Gorithm to Solve Non-positive-definite Fi- Nite Element Equations in Electromagnetic Problems
نویسندگان
چکیده
A new method called Expanded Cholesky Method (ECM) is proposed in this paper. The method can be used to decompose sparse symmetric non-positive-definite finite element (FEM) matrices. There are some advantages of the ECM, such as low storage, simplicity and easy parallelization. Based on the method, multifrontal (MF) algorithm is applied in non-positive-definite FEM computation. Numerical results show that the hybrid ECM/MF algorithm is stable and effective. In comparison with Generalized Minimal Residual Method (GMRES) in FEM electromagnetic computation, hybrid ECM/MF technology has distinct advantages in precision. The proposed method can be used to calculate a class of non-positivedefinite electromagnetic problems.
منابع مشابه
Efficient Multigrid Solvers for the Stokes Equations using Finite Elements
Mixed nite element application to solve the Stokes equations using di erent smoothers by mean of multigrid technique has been investigated. Multigrid technique's e ciency for positive de nite linear systems is proved but it is much more challenging in the case of saddle-point problems. The e ciency of a multigrid method is highly dependent on the smoother and the coarse grid solver. A two-dimen...
متن کاملA Numerical Method for Solving the Elliptic Interface Problems with Multi-domains and Triple Junction Points
Elliptic interface problems with multi-domains and triple junction points have wide applications in engineering and science. However, the corner singularity makes it a challenging problem for most existing methods. An accurate and efficient method is desired. In this paper, an efficient non-traditional finite element method with non-body-fitting grids is proposed to solve the elliptic interface...
متن کاملA new positive definite semi-discrete mixed finite element solution for parabolic equations
In this paper, a positive definite semi-discrete mixed finite element method was presented for two-dimensional parabolic equations. In the new positive definite systems, the gradient equation and flux equations were separated from their scalar unknown equations. Also, the existence and uniqueness of the semi-discrete mixed finite element solutions were proven. Error estimates were also obtaine...
متن کاملIMF: An Incomplete Multifrontal LU-Factorization for Element-Structured Sparse Linear Systems
We propose an incomplete multifrontal LU-factorization (IMF) preconditioner that extends supernodal multifrontal methods to incomplete factorizations. It can be used as a preconditioner in a Krylov-subspace method to solve large-scale sparse linear systems with an element structure; e.g., those arising from a finite element discretization of a partial differential equation. The fact that the el...
متن کاملAn Efficient Strain Based Cylindrical Shell Finite Element
The need for compatibility between degrees of freedom of various elements is a major problem encountered in practice during the modeling of complex structures; the problem is generally solved by an additional rotational degree of freedom [1-3]. This present paper investigates possible improvements to the performances of strain based cylindrical shell finite element [4] by introducing an additio...
متن کامل